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Abstract. The tremendous experimental progress in cosmic microwave background (CMB) temperature and 
polarization anisotropy studies over the last few years has helped establish a standard paradigm for cosmology at 
intermediate epochs and has simultaneously raised questions regarding the physical processes at the two opposite 
ends of time. What is the physics behind the source of structure in the universe and the dark energy that is 
currently accelerating its expansion? We review the acoustic phenomenology that forms the cornerstone of the 
standard cosmological model and discuss internal consistency relations which lend credence to its interpretation. 
We discuss future milestones in the study of CMB anisotropy that have implications for inflationary and dark 
energy models. These include signatures of gravitational waves and gravitational lensing in the polarization flelds 
and secondary temperature anisotropy from the transit of CMB photons across the large-scale structure of the 
universe. 



1. INTRODUCTION 

The pace of discovery in the field of Cosmic Microwave Background (CMB) anisotropy has been accelerating 
over the last few years. With it, the basic elements of the cosmological model have been falling into place: 
the nature of the initial seed fluctuations that through gravitational instability generated all of the structure 
in the universe, and the mixture of matter-energy constituents that drives its expansion. Since even reviews 
written a scant year ago (e.g. ||, D) seem out of date today, we will focus mainly on the fundamental 
physical elements that enter into cosmic microwave background observables. 

In 1992, the COBE DMR experiment reported the first detection of cosmological anisotropy in the 
temperature of the CMB The 10~^ variations in temperature detected on scales larger than the 7° 
resolution provided strong support for the gravitational instability paradigm. These variations represent the 
direct imprint of initial gravitational potential perturbations through their redshifting effect on the CMB 
photons, called the Sachs- Wolfe effect and are of the right amplitude to explain the large-scale structure 
of the universe. From 1992-1998, a host of experiments (see references in [||) detected a rise and fall in 
the level of anisotropy from degree scales to arcminute scales. In the following two years, the Toco |], 
Boomerang and Maxima ||] experiments measured a clearly defined first peak in the spectrum and 
provided empirical evidence that the small scale anisotropy is dominated by coherent acoustic phenomena, 
long predicted to exist jl^, |ll|. This first peak measurement has subsequently been confirmed by several 
groups, with the best measurement to date from the Archeops experiment | ]l2| . These measurements firmly 
establish initial gravitational potential or curvature fluctuations as the primary source of structure in the 
universe (c.f. cosmological defect models |lj, |l^), provide clear evidence that the universe is close to 
spatially flat and, in conjunction with other cosmological measurements of the dark matter, an indicate that 
there is a component of missing or dark energy. These findings provide support for the inference from distant 
supernovae that the expansion is currently accelerating under the influence of the dark energy ]l6| , . 

In 2001, the DASI and Boomerang ||l^, ^ experiments announced a detection of secondary acoustic peaks 
in the spectrum which was subsequently confirmed by the VSA experiment . These experiments provided 
a precise measurement of the baryon-photon ratio that is in excellent agreement with inferences from big 
bang nucleosynthesis (e.g. pH), lending further confidence in the underlying cosmological model. Moreover 
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they provided the first direct evidence that dark, non-baryonic, matter exists in the early universe. In the 
past year, two more predicted phenomena have been discovered: the dissipation of the acoustic peaks at 
small scales by the CBI experiment and, following a series of increasingly tighter upper limits 
the polarization of the anisotropy by the DASI experiment [ p5| . These observed phenomena provide the best 
internal evidence that the physical assumptions underlying the interpretation of the peaks are justified. 

The past few years have been an era of milestones for experiments and millstones for theoretical speculation. 
The standard cosmological paradigm of structure formation by gravitational instability of cold dark matter 
in a nearly homogeneous and isotropic universe has been so thoroughly tested that it is extremely difficult 
to find viable contenders that differ in any fundamental way. The acoustic peaks are furthermore strikingly 
consistent with two fundamental predictions of simple infiationary models of the early universe: a fiat spatial 
geometry and a nearly scale-invariant spectrum of initial curvature fluctuations. Alternatives to inflation, 
must now be nearly indistinguishable phenomenologically to be consistent with the data (e.g. [^). 

The future may see a reversal of this trend as CMB experiments and theory focus on mysteries at the two 
opposite ends of time. The detailed physics of inflation (or contending theories) and the dark energy remains 
unexplained and largely unexplored. Here the increase in precision, and perhaps more importantly accuracy, 
expected from future experiments will allow studies of the associated phenomena: the statistical properties of 
the anisotropy, features in the initial curvature spectrum, subtle effects in the polarization, and the secondary 
anisotropy generated by structure in the universe. In fact secondary anisotropy from unresolved clusters of 
galaxies may have been detected recently by the CBI and BIMA experiments. 

Key past and future milestones are summarized in Table 1. In this review we shall try to place this 
somewhat bewildering array of phenomena in its cosmological context. We begin with a description of CMB 
observables in §^ proceed through the physical basis of their generation in to the phenomenology of 
the acoustic peaks and secondary anisotropy in §§j|. Finally in §^ we outline future milestones of CMB 
temperature and polarization anisotropy. 



2. CMB OBSERVABLES 

The fundamental observable in the CMB is the intensity of radiation per unit frequency per polarization at 
each point in the sky. The polarization state of the radiation in a direction of sky denoted n is described 
by the intensity matrix (^Ei{n)E* {nj'j where E is the electric field vector and the brackets denote time 
averaging. 

The radiation is measured to be a blackbody of T = 2.728 ±0.004K (95% CL) [P9[ with fractional variations 
across the sky at the 10"^ level and fractional polarization at the 10"® level [^5[. It is then convenient to 
describe the observables by a temperature fluctuation matrix decomposed in the Pauli basis (e.g. 0|) 



P = C(E(n)Et(n)) 

= e(n)I + Q(n) (73 + U{h) cTi + V{n) , 



(1) 




FIGURE 1. Power spectra data plotted as the rms contribution per logarithmic interval + l)Ci/2n]^^^ with error boxes 
representing Icr error bars and approximate multipole bandwidth. Overplotted is a scale-invariant, flat cosmological model with 
f^m = 1/3, Ha = 2/3, h = 0.7, f2i,/i^ = 0.02, reionization redshift Zi = 7, and an inflationary energy scale of Ei = 2.2 X lO^'' GeV. 



where we have chosen the constant of proportionahty so that the Stokes parameters {<d,Q,U,V) are dimension- 
less, e.g. 6(n) = Ar(n)/T averaged over polarization states. Note that the circular polarization V is absent 
cosmologically and under a counterclockwise rotation of the coordinate axes by -0, Q±iU ^ e^'^''^ {Q ± ill) . 
The temperature and polarization fields are decomposed as ^ 

e,™ = J dnY;^{h)e{h), 

E,^ ± iBi„, = -J dn±2 YL (n) ± if/(n)] , (2) 

in terms of the complete and orthogonal set of spin harmonic functions, sY/m, which are eigenfunctions of the 
Laplace operator on a rank s tensor p^ , |35[ |. Yim =^ o^im is the ordinary spherical harmonic. For small sections 
of sky, the spin-harmonic expansion becomes a Fourier expansion with Yim — > e'' " and ±2Yim — *■ — e='=^*'^'e'' ", 
where (j>i is the azimuthal angle of the Fouler wavevector 1. Note that the E and B modes are then simply 
the Q and U states in the coordinate system defined by 1 p6| , i.e. -B-modes have a polarization orientation at 
45° to the wavevector. This decomposition is analogous to the vector field case where the divergence-free and 
curl-free portions are distinguished by the orientation of the velocity vector with respect to the wavevector. 
Represented as E{n) = J^^imYimifi) and B{n) = J2BimYim{f>), they describe scalar and pseudoscalar fields 
on the sky and hence are also distinguished by parity. 

The primary observable is the two point correlation between fields X,X' e {ld,E,B} 

{^lra^l'in')=^U'^rnm'C^^ , (3) 

and are described by power spectra C/ as long as the fields are statistically isotropic. If parity is also conserved 
then B has no correlation with O or E. If in addition the fiuctuations are Gaussian distributed, the power 
spectra contain all of the statistical information about the fields. The measurements of these power spectra 
to date are shown in Fig. H^. 
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FIGURE 2. Geometric projection of sources at the last scattering surface into temperature and polarization anisotropy 
today. Left: the main sources of anisotropy are the effective local temperature which has a monopole structure (yoo)i the Dopplcr 
effect which has a dipole structure (Vio) and the quadrupole source of polarization {Y20) decomposed into plane waves. Right: 
the projection of the sources involves a coupling of the plane wave angular dependence or momenta with that of the source. A 
plane wave temperature inhomogeneity produces a sharp peak in the anisotropy at I ft: kDt{= 100 here) reflecting the structure 
of j;(fcD») shown; the Doppler effect does not produce peaks since the coupling involves i±l in the combination jj(fcD«). The 
functions here are weighted by 2i + 1 to account for the number of modes. 



3. PHOTON-BARYON PLASMA 
3.1. Cosmological Background 

We begin with a brief review of the cosmological background and the parameters that govern it. The 
expansion of the universe is described by the scale factor a(i), set to unity today, or equivalently the redshift 
z{t) = a(t)~^ — 1, and by the current expansion rate, the Hubble constant _ffo = lOOft. km sec~^ Mpc^^. The 
universe is flat (no spatial curvature) if the total density is equal to the critical density, pc = l.SSh? x 10~'^^g 
cm^''; it is open (negative curvature) if the density is less than this and closed (positive curvature) if greater. 
The mean densities of different components of the universe control a{t) and are typically expressed today in 
units of the critical density $1^, with an evolution pi cx a"'^'-^^'"'-*, where Wi =pi/ pi with pi is the pressure. In 
particular the CMB photons have — l/'i and so p^ cx a^^ or Toe a"^ . The quantities ^lih? are proportional 
to the physical density of the species today. We will be interested in the baryonic component Q.^h? and the 
total non-relativistic matter Q.„ih? , where the difference is made up of cold dark matter. We also allow for a 
dark energy component fie, with an equation of state parameter We which for illustrative purposes we often 
take to be a cosmological constant (f^A, WK = ^1) as consistent with indications from distant supernovae 
| p6| , p7[ |. We represent the total density as fitot = X^i^i so that deviations from unity represent a spatially 
curved universe. 



3.2. Sources of Anisotropy 

As the CMB cools due to the expansion, the photons eventually no longer have sufficient energy to overcome 
the binding energy of hydrogen B = 13.6eV to keep the medium ionized. After this epoch, CMB photons 
propagate essentially unimpeded to the observer. This transition, called recombination, occurs relatively 
rapidly around a* « lO"-^ at an energy scale of 1/3 eV as can be seen by the equilibrium Saha equation 

x1 rieUp 1 / TOgT 

1 — Xe nnrib rib \ 2tt 

V 0.02 \T 




where Xg is the ionization fraction of hydrogen and we have neglected a small contribution from helium. 
The low T/B of the transition is mainly due to the low baryon-photon ratio of the universe. It is sufficient 
that photons in the tail of the blackbody distribution be energetic enough to ionize hydrogen. In reality, 
recombination is slightly less rapid than implied by the Saha equation since processes which account for net 
recombination are not sufhciently rapid to maintain equilibrium toward the end of recombination Q (see 
p9[ for refinements). Nonetheless the transition is sufficiently rapid that the photons in the CMB originate 
from a sharply defined surface or photosphere called the last scattering surface. 

The properties of the plasma on this surface directly translate into the observable primary temperature 
and polarization anisotropy in the CMB. These properties are themselves governed by Thomson scattering 
of photons off of free electrons, which has a differential cross section of 

|^ = ^|E'.E|V., (5) 

where ctt — Stto^ /Sm,, is the Thomson cross section, E' and E denote the incoming and outgoing directions 
of the electric field or polarization vector. Summed over polarization and incoming angle, the comoving mean 
free path of a photon in the electrons is 



(^)'mpc, (6) 



Ueara \ 0.02 J V 10 

which is relatively small by cosmological standards. Overdots here and below represent derivatives with 
respect to conformal time. Electrons are coupled to baryons by Coulomb interactions and so on scales 
larger than r~^, the photon-baryon plasma can be considered a nearly perfect fluid ]lO| . In particular, rapid 
scattering keeps the photons isotropic in the baryon rest frame and so sets its dipole moment equal to the 
baryon velocity. 

The tight coupling approximation implies that the CMB is described by relatively few quantities. There 
is the local temperature, or monopolc, in space Ooo(x), the dipole and a residual but important 

quadrupole moment 62m (x) from the breakdown of tight coupling. In addition, the photons experience 
a gravitational redshift from gravitational potentials and gravitational waves. Modern CMB codes ^ 
track the evolution of these sources ^ and construct the present day CMB temperature and polarization 
anisotropy spectrum by geometric projection. 

As a simple example, let us consider the spatial temperature perturbation 0(x) = '\/47r0oo(x), where the 
normalization convention reflects that of Yim, on a last scattering surface considered to be infinitely sharp 
at a comoving distance (see Fig. |^). In a flat universe, the spatial field may be represented by its Fourier 
harmonics 



e(n) = J dDe{x)6{D~D^ 



'^''^ e(k)e^''-^*". (7) 



(27r 



l3 



In an spatially curved universe, the plane waves must be replaced by the eigenfunctions of the Laplace 
operator to account for a change in the relationship between distance and angle that we discuss below. Also 
to achieve high precision in the predictions, the delta function must be replaced by the visibility function 
Te""^, the probability of a photon last scattering in a distance interval dD to account for the finite duration 
of recombination. 

Expanding the plane waves in spherical harmonics, we obtain 

/d^k 
^^e(k)Wj,(fci?,)rz„(k). (8) 

The power spectrum of the anisotropy contributions from the local temperature is then governed by that of 
the spatial temperature field at recombination 



(e*(k)e(k')) = (2^)35(k-k')PT(fc), 



(9) 



such that 



cr-^^JfjUkD.)'^. (10) 

Note that the quantity fc^Pr/STr^ is the contribution per log interval to the variance of the temperature 
(T^ = J d'^ k / {271)^ Pt ■ For a slowly-varying log power, the integral in Eqn. ( |lO| ) can be performed analytically 

and so it is convenient to represent the anisotropy by the rms temperature contribution per log interval 
At- It is also the contribution to the anisotropy variance per log interval in / for 1^1. We also define 
the analogous quantities to describe the polarization fields (Ap) and the temperature-polarization cross 
correlations {l{l + l)Cf^ /2Tr). 

This behavior of the local temperature is representative of the other sources with a few important 
differences. In general, power from low order multipole moments I = 0,1,2 in the CMB project into higher 
order anisotropy today by simple geometry The local temperature is an / = source. Note that the 
spherical Bessel function ji{x) peaks at / ~ cc (see Fig. || right) and so the translation in Eqn. (^) from power 
in real space to angular space is nearly one-to-one. Thus in Eqn. (|l^) one takes « I with features in one 
preserved in the other. This corresponds to an edge-on (n _L k) projection of the wavelength of the fluctuation 
onto angle: 9 = X/D^. 

Other sources of anisotropy have angular dependence on the last scattering surface which changes the 
relationship between I and k. The dipole term provides an instructive example. It represents Doppler shifts 
in the temperature due to the motion of the photon-baryon plasma along the line of sight. For a curl free flow 
as expected cosmologically, there is no contribution to the Doppler effect for n _L k and the peak structure 
is washed out. Mathematically, the Doppler effect contributes a dipole term fi - k oc Yi„i which when coupled 
to the plane wave angular momentum Yi„i leads to a recoupling of ji±i to j; (see Fig. Unlike ji, j[ has 
no peak structure leading to a one-to-many projection of k onto I. Geometrically, well-defined features such 
as peaks in the anisotropy spectrum cannot come from the Doppler effect pif . All sources may be described 
in this framework of the addition of local angular momentum with plane wave orbital angular momentum 
with the net effect of modifying the mapping between k and I space . 

The local quadrupole plays a special role in that it is the fundamental source of linear polarization in 
the CMB as can be seen from the Thomson differential cross section of Eqn. (^. Heuristically, incoming 
radiation shakes an electron in the direction of its electric field vector E' causing it to radiate with an outgoing 
polarization parallel to that direction. However since the outgoing polarization E must be orthogonal to the 
outgoing direction, incoming radiation that is polarized parallel to the outgoing direction cannot scatter, 
leaving only one polarization state. If the intensity were completely isotropic the missing polarization state 
is supplied by radiation incoming from the direction orthogonal to the original one. Only a quadrupole 
temperature anisotropy in the radiation generates a net linear polarization from Thomson scattering. In 
particular, the polarization is oriented along the cold axis of the quadrupole moment in the transverse plane. 
We shall see that the breakdown of tight coupling during recombination allows quadrupoles to be generated 
out of gradients in the fluid velocity and from gravitational redshifts associated with gravitational waves. 
For the former, since the quadrupole axis is aligned with the wavevector, the polarization generated is a pure 
i?-mode (see Fig. ||). Gravitational waves generate both E and B modes 

To understand the phenomenology of the anisotropy in the CMB temperature and polarization fields it is 
sufficient to understand the dynamics governing the local sources: the monopole, dipole and quadrupole of 
the CMB temperature field and their gravitational or metric fluctuation analogues. 



3.3. Fluid Dynamics 

The evolution of the anisotropy sources is governed by simple fluid dynamics. From the previous section, 
the critical variables are the monopole or temperature fluctuation 8(k) = V47rOoo(k), the dipole or bulk 
velocity ^^(k) = — i-\/47r/30io(k), and the quadrupole or anisotropic stress 7r^(k) — —( 12/5) 1/471/5 G2o(k). 



k=l/D. 



(11) 



We will hereafter drop the argument k with the understanding that the Fourier representation is always 
assumed. For convenience we have chosen the coordinate system so that z || k so that the plane waves are 
azimuthally symmetric and stimulate only the m — mode (see Fig. |^). Likewise we have suppressed the 
vector dependence of the bulk velocity by the same assumption = — iu^k. The analogous quantities for the 
baryons are the density per turbation Sf, and bulk velocity Vb. For gravity, we choose a conformal Newtonian 
representation (see e.g. P6|| ) where the gravitational potential perturbations are defined by the Newtonian 
potential '5 (time-time metric fluctuation) and the curvature fluctuation $ (space-space metric fluctuation 

Covariant conservation of energy and momentum requires that the photons and baryons satisfy separate 
continuity equations 



e^-^v^-^, 6b^-kvb-3^, (12) 



and coupled Euler equations 



=: fc(e + ^') - ^TTy ~ f{v^ - Vh) , 

Vb^-^Vb + k'i' + T{vy-Vb)/R, (13) 

where R = {pb + pb) / {p-y + Pj) ~ ipb/'^p^ is the photon-baryon momentum density ratio. We have neglected 
a small correction to the anisotropic stress term in a curved universe. 

The continuity equations represent particle number conservation. For the baryons, pb oc n^. For the photons, 

1 /3 

T cx n^l , which explains the 1/3 in the velocity divergence term. The <f> terms come from the fact that $ 
is a perturbation to the scale factor and so they are the perturbative analogues of the cosmological redshift 
and density dilution from the expansion. The Euler equations have similar interpretations. The expansion 
makes particle momenta decay as a^^ . The cosmological redshift of T accounts for this effect in the photons. 
For the baryons, it becomes the expansion drag on Vb {a /a term). Potential gradients k^ generate potential 
flow. For the photons, stress gradients in the fluid, both isotropic {kSp^/{pj + pj) = kQ) and anisotropic 
{kiTy) counter infall. Thomson scattering exchanges momentum between the two fluids (f terms). 

For scales much larger than the mean free path f~^, the Euler equation may be expanded to leading order 
in fc/r, such that the photons are isotropic in the baryon rest frame Vj = Vb and so the joint Euler equation 
becomes 

^[{1 + R)v.,]^k[e + {1 + R)^]. (14) 
ai] 

Combining this with the continuity equation leads to the oscillator equation (e.g. Q]) 

±[ii+R)e] + ^e = -'^ii + R)^--^[ii + Rm, (15) 

and a small residual anisotropic stress or quadrupole that tracks the evolution of the fluid velocity , 

7r^ = — -w^. (16) 

iO T 

This dependence reflects the fact that a local quadrupole can arise from a gradient in the velocity fleld, for 
example as photons from two hot crests of a plane wave fluctuation meet at the trough in between (see Fig. ^ 
left), but is suppressed by scattering. 

Equation (|l^) is the fundamental relation for acoustic oscillations. The change in the momentum of the 
photon-baryon fluid is determined by a competition between the pressure restoring and gravitational driving 
forces which causes the system to oscillate around its equilibrium. Note that the frequency of the oscillation 



where Cg is the sound speed of the fluid. 



Equation (16) is the fundamental relation for acoustic polarization. Its quadrupole source tracks the 
motion of the fluid. Polarization is generally small, at most ~10% of the anisotropy itself, since it requires 
scattering for its generation but its quadrupole source is suppressed if the scattering is too rapid. Given the 
initial conditions and gravitational potentials, these equations predict the phenomenology of the acoustic 
oscillations in the temperature and polarization fields. 



3.4. Initial Conditions 

The simplest inflationary models make a set of definite predictions for the initial conditions of the acoustic 
oscillations and hence their successful observation provides strong support for the inflationary paradigm. 
Quantum fluctuations in the scalar fleld that drives inflation imprints a nearly spectrum of Gaussian curvature 
(potential) fluctuations fc^P$/27r^ c>c k"~^ where n w 1 |4^, |5^ on a spatially flat background metric. 
Gravitational infall into these initial potentials eventually generates all of the structure in the universe. 
Quantum fluctuations in the gravitational wave degrees of freedom also produce a nearly scale-invariant 
spectrum of fluctuations whose power depends on Ef where is the energy scale of inflation ^ . 

A Newtonian gravitational potential w — $ necessarily imparts an initial temperature perturbation 
since ^' represents a spatially varying time-time perturbation to the metric away from coordinates where the 
temperature is homogeneous. The perturbation is equivalent to a change in the scale factor since 

acxt3(i+p/p) ^ (18) 
which then produces a change in the temperature perturbation from the cosmological redshift T oc a^^ of 

or — 5'/2 in the radiation dominated era |^, We call 6 + ^ the effective temperature since it also accounts 
for the redshift a photon experiences when climbing out of a potential well, also known as the Sachs- Wolfe 
effect In the matter dominated epoch, O -I- = ^f/S. 

There are three important aspects of these results. First, inflation sets the temporal phase of all wavemodes 
by starting them all at the initial epoch. We shall see that this predicts a coherent set of peaks in the 
CMB spectrum with a definite phase. Cosmological defect models predict a more random distribution of 
acoustic phases which produces incoherent acoustic phenomena |l^ . Defects can be now be ruled out 

as a primary mechanism for structure formation in the universe. More generally, without inflation or some 
other modiflcation to the matter-radiation dominated universe, curvature perturbations cannot be generated 
outside the apparent horizon and so build up only by the causal motion of matter. This generally entails at 
least a delay in temporal phase of the oscillations which is not observed. 

Secondly, since the power spectrum of the effective temperature is directly related to scale-invariant 
curvature fluctuations from inflation, Eqn. ( pT[ ) implies the acoustic oscillations should also be approximately 
scale invariant in amplitude in the tight coupling regime. Observations are in excellent agreement with this 
fundamental prediction with tight constraints on the index n = 0.941q |!)4 (|55 , see also 56, 5^, Q). We will 



therefore base the discussion of the acoustic phenomenology on models with nearly scale-invariant initial 
curvature fluctuations as predicted by inflation. 

Finally, the scale-invariant gravitational wave background leads to a quadrupolar distortion in the CMB 
temperature field just like its effect on a ring of test masses. Because the quadrupole axes lie in the plane 
transverse to the wavevector, this quadrupole anisotropy leads to i?-mode polarization as the tight coupling 
approximation breaks down. A measurement of i?-modes from gravitational waves would determine the 
energy scale of inflation Ei but its strong scaling with Ei implies that it can only open a relatively small 
window between a few 10^^ - few 10^^ GeV for possible detection l59|. 
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FIGURE 3 

model of Eqn. (|23D 

same time r] = -q 



bmporal evolution of the effective temperature as a function of wavenumber predicted by the simplified 
Left: amplitude; right: rms fluctuation. All modes begin with the same phase and are frozen in at the 
The result is a harmonic series of extrema with a fundamental scale of = 7r/s». Baryons, here shown with 



an i? = 0.6{?;/»7,)^, displace the zero point of the oscillation and enhance the compressional troughs. In the rms, this leads to 
enhanced odd numbered peaks. 



4. ACOUSTIC PHENOMENOLOGY 
4.1. Acoustic Scale 

It is instructive to first consider a simplified model where the universe is always matter-dominated in its 
expansion and the dynamical effects of the baryons are negligible. This amounts to holding $ and constant 
in the oscillator equation (^5|) and setting i? = 0. The solution for an initial curvature fluctuation is then 
simple 

[e + «'](/s,r/) -^'(fc,0)cos(fcs), i;^(fc,?7) = ^*(fc,0)sin(fcs), (20) 

where s — Csdrj' is the distance sound can travel by 7] or the sound horizon. 

In the limit of scales large compared with the sound horizon fcs <C 1, the perturbation is frozen into 
its initial conditions. This is the gist of the statement that the large-scale anisotropy measured by COBE 
directly measure the initial conditions. On small scales, the amplitude of the Fourier modes exhibits temporal 
oscillations corresponding to compression and rarefaction of the plasma inside gravitational potential wells. 
All modes are frozen in at recombination 77* (see Fig. ^ left). Modes that are caught at either maxima or 
minima of their oscillation at recombination correspond to peaks in the rms fluctuation or power (see Fig. ^ 
right). Given the solution in Eqn. (|2^), these modes are harmonics fc„ = nfc^ of a fundamental acoustic 
scale kA = ti'/s*, where n is an integer (see Figure ^). From Eqn. (|ll|), these become angular peaks in 
the anisotropy power with a characteristic angular scale of I a — fcyiZ?,. The Doppler effect from the line- 
of-sight velocity has an rms v~^l\fZ and so contributes an equal amplitude fluctuation. It is 7r/2 out of 
phase with the temperature. This is because extrema represent turning points in the oscillation where the 
velocity vanishes. Nonetheless, because of the geometry of the projection (see Fig. these oscillations do 
not contribute peaks in the angular power spectrum. We shall see that they do however predict the peaks 
in the polarization spectrum through Eqn. ([16[). 

Under the flat, matter-dominated assumption the horizon distance scales as ry oc a}^!"^ . Thus w Z?* /a/3000 
and so Ia ~ 200. In a spatially curved universe, the distance used in the conversion from fc to / no longer 
equals the comoving distance -D* 7^ 77(0 = 1) — ry*. Consider first a closed universe with radius of curvature 
i? = 7Jjj'^|f2tot — Suppressing one spatial coordinate yields a 2-sphere geometry with the observer 

situated at the pole (see Fig. ^ left). Light travels on lines of longitude. A physical scale A at fixed latitude 
given by the polar angle Q subtends an angle a = X/Rsind. For a <C 1, a Euclidean analysis would infer a 
distance D = RsinO, even though the coordinate distance along the arc is d = OR; thus 



D = Rsin{d/R) 



(21) 
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FIGURE 4. The acoustic scale. Left: geometrical effect in a closed universe. Objects in a closed universe are further than 
they appear {D < d). Consequently the measured angular scale of the peaks is extremely sensitive to the curvature, once the 
physical scale has been set by the acoustic physics. Right: translation of a constraint on the peak scale 300 <Ia < 308 fed] (68% 
CL) onto the cosmological parameters Qrn and Q\ with alLqther parameters fixed versus constraints with uncertainties in the 
physical matter density Qmh^ from CMB determinations |55| folded in. Uncertainties in this plane are no longer dominated by 
peak location measurement errors. 



For open universes, simply replace sin with sinh. The peak locations In = nkAD^ making them extremely 
sensitive to the spatial curvature JTl] , |6C| |; measurements constrain the geometry to be close to flat, consistent 
with a radius of curvature that is larger than the observable universe. Since a flat universe is at the critcal 
density, local measurements of a low matter density indicate that there is a missing or "dark" component of 
the energy density, in good agreement with indications of an accelerating expansion from distant supernovae. 

These simple scalings must be modified for the fact that the universe is not fully matter dominated at 
recombination, the fluid is not purely photon-dominated, and the expansion is dark energy dominated today. 
The main effect is from the radiation density which changes the expansion rate and so shifts the acoustic 
scale to higher multipoles. For a flat universe with flmh'^ ~ 0.15, Ia ~ 300 which is a substantial shift. The 
baryons lower the sound speed and so have a similar but smaller effect. The dark energy also provides a 
smaller effect through a decrease in the distance to recombination and hence a lowering of I a, with an 
increase in fl^ or w. The sensitivity of the acoustic scale to cosmological parameters is approximately (see 
also Fig. H) 

^ -^tot m _Q Yj_^ _Q Y\Aw^ + ^Sf! \ , (22) 



around a model of fitot — 1, ^mh^ — 0-15, J7e — 0.65, We — —1 and ^i,h^ — 0.02 A joint analysis of the 
data yields I a = 304 ±4 and with these errors, the uncertainty in the cosmological interpretation in the 
dark energy and curvature domain is already dominated by uncertainty in Q,mh? not peak measurement error 
(see Fig. ^ right). Fortunately as we shall see, this parameter can be internally determined from the CMB 
and serves as an example where future increased precision will have important implications for a quantity of 
fundamental interest, the dark energy. 

Finally, there is a separate effect on the first peak. We shall see that the decay of the gravitational potential 
in a radiation dominated universe generates anisotropy from gravitational redshifts after last scattering, filling 
in the rise to the first peak and shifting its location downwards off of the harmonic series to li ~ SI a/4:, placing 
li « 220 in agreement with the observed location p^ . 



4.2. Baryon loading 



Baryons add to the mass of the photon-baryon plasma without adding to the pressure. An examination 
of Eqn. (15) shows that their effect comes solely through the baryon-photon momentum density ratio 
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FIGURE 5. Sensitivity of the temperature power spectrum to four fundamental cosmological parameters: the energy 
density of the dark energy today Qe in units of the critical density, the equation of state parameter of the dark energy We, the 
physical baryon density Hj/i^ and the physical matter density Q,mh^ ■ AH are varied around a fiducial fiat model of f2e = 0.65, 
Wf, = -1, n^/i^ = 0.02 and Vlmh? = 0.15 with n = 1. 



R fn 0.6(ilhft-^/0.02)(a/10^'^). It is instructive to look at the solution to the oscillator equation ( |2(]| ) in the 
approximation that R is constant. Again under the assumption of a matter-dominated universe, the solution 
is that of a simple harmonic oscillator in a constant gravitational field but with an increased mass term 

[6 + *] (A:,?/*) = [e + (l + i?)*](fc,0)cos(A:s) 

= [l + 3i?]i*(0)cos(fcs)-i?^'(0), (23) 
o 

where the sound speed entering into the sound horizon calculation is reduced since cj. = 1/3(1 + R). Aside 
from this reduction, baryons have two distinguishing effects: they enhance the amplitude of the oscillations by 
l-|-3i? fractionally and shift the zero point of the oscillation by —R'^. The latter modulates the amplitude 
of neighboring peaks: odd numbered peaks will be enhanced over the zero-baryon case by 1 + 6i?; even 
numbered peaks remain the same (see Fig. right). Physically, the baryon mass enhances compression 
inside gravitational potential wells. 

These qualitative results remain true in the presence of the real time-variable R. Measurement of these 



baryonic signatures currently limits fibft,^ ~ 0.022lg'QQ2 |55 . This value is strikingly consistent with inferences 
of the baryon density from big bang nucleosynthesis. Consequently there have been no significant changes in 
the baryon-photon ratio from an energy scale of an MeV to an eV in the expansion history of the universe. 



4.3. Dark Matter and Radiation 



The matter-to-radiation ratio scales as Pm/ Pr ~ 3.6(il„i^^/0.15)(a/10^'^) and so the universe is only barely 
matter-dominated at last scattering. Moreover fluctuations corresponding to the higher peaks entered the 
sound horizon at an earlier time, during radiation domination. As we have seen, including the radiation 
changes the expansion rate of the universe and hence the physical scale of the sound horizon at recombi- 
nation. Fortunately, radiation has the more unique effect of driving the acoustic oscillations by making the 
gravitational force evolve with time [Q . 

The exact evolution of the potentials is determined by the relativistic Poisson equation. Its qualitative 
content is clear: since the background density is decreasing with time, the density fluctuations must grow 
unimpeded by pressure to maintain constant potentials. In particular, in the radiation dominated era 
once pressure begins to fight gravity at the first compressional extrema of the oscillation, the Newtonian 
gravitational potential and spatial curvature must decay. 

This decay actually drives the oscillations: it is timed to leave the fluid maximally compressed with no 
gravitational potential to fight as it turns around. The net effect is doubled since the redshifting from the 
spatial metric fluctuation $ also goes away at the same time. When the universe becomes matter dominated 
the gravitational potential is no longer determined by photon-baryon density perturbations but by the 
pressureless dark matter. Therefore, the amplitudes of the acoustic peaks increase as the matter-to-radiation 



ratio decreases |62, The net result is that across the horizon scale at matter-radiation equality the 
acoustic amplitude increases by a factor of 

^^ = 5 (24) 

3 ^ 

for a pure photon and dark matter universe, and a factor of 4 when including the effect of neutrinos and 
baryons |Q (see also Fig. ||). By eliminating gravitational potentials, radiation also eliminates the alternating 
peak heights from baryon loading (see Fig. |^). The observed high third peak (see Fig. ^ is a good indication 
that matter dominates the energy density at recombination. Finally, the effect of the decaying potential 
after recombination leads to so-called integrated Sachs Wolfe contributions to the temperature fluctuations 
through the continuity equation (|l^) and shifts the flrst acoustic peak downwards off of the acoustic series 

il- 

Observations of these matter-radiation phenomena (see Fig. ||) currently constrain the total matter density 
to be ^mh'^ = 0.15 ±0.03 |^ and are crucial in internally resolving ambiguities in the interpretation of the 
peak scale. Since this number greatly exceeds the measured baryon density, these observations are also the 
first empirical indication that the dark matter exists at high redshift. 



4.4. Damping 

The photon-baryon fluid has slight imperfections corresponding to shear viscosity ir-y and heat conduction 
(vb — Vj) in the Euler equation (|l^) ^4, |6^. These imperfections damp acoustic oscillations. These are 
both associated with the diffusion of the photons which is especially important during recombination 
when the mean free path of the photons rises dramatically. As we have seen the viscosity tTj ~ kv^/i. 
With the continuity equation kv^ « —30, this leads to a damping term in the oscillator equation. The 
heat conduction term can be shown to have a similar effect by expanding the Euler equations in k/f. 
An examination of the resulting damped oscillator equation shows that diffusion suppresses the oscillation 
amplitude by a factor of order The damping scale fed is thus of order \/f /rj, corresponding to the 

geometric mean of the horizon and the mean free path. Detailed numerical integration of the equations 
of motion are required to track the rapid growth of the mean free path and damping length through 
recombination itself. These calculations show that the damping scale is of order fcdS* « 10 leading to a 
substantial suppression of the oscillations beyond the third peak (see Fig. H). Observations of the damping 
phenomena provide a check on the fundamental assumptions underlying the interpretation of the acoustic 
peaks. With Q.i,h^ and flmh^ and the acoustic scale Ia flxed by the first few peaks, the damping scale is 
uniquely predicted given the atomic physics of recombination. Any change in recombination, for example due 
to a variation in the fine structure constant, will be revealed as a discrepancy in the predictions. Likewise, 
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FIGURE 6. Phenomenological contributions to the anisotropy form. Radiation drives the amplitude of the acoustic 
oscillations up across the multipole scale subtended by the horizon aLmatter radiation equality /oq and down across the 
damping scale Ij^. Functional forms for these envelopes are taken from [pM]. The ISW effect pushes l\ <Ia due to the radiation 
and raises the large angle anisotropy due to the dark energy. The polarization scales as Ijlj:, for I Kip due to the tight-coupling 
suppression of its quadrupole source and is also damped below I u . 



any misinterpretation of I a due to say initial conditions that were not purely curvature fluctuations resulting 
in a phase shift in the peaks would also show up as a discrepancy. The data in the damping tail to date are 
beautifully consistent with the predictions of the model (see Fig. 



4.5. Polarization 

The dissipation of the acoustic oscillations leaves a signature in the polarization of the CMB in its wake. 
Recall that the observed polarization is a projection of the local quadrupole temperature anisotropy at last 
scattering. The fact that the polarization source is the quadrupole explains the shape and height of the 
polarization spectra in Fig. ^ Since the quadrupole is of order kv^/r ^ (k/k^i){kd'q^)~'^v^ (see Eqn. 16) 



the polarization spectrum rises as I /Id to peak at the damping scale with an amplitude of about 10% 
of the temperature fluctuations before falling due to the elimination of the acoustic source itself due to 
damping. Since is out of phase with the temperature, the polarization peaks are also out of phase with the 
temperature peaks. Furthermore, the phase relation also tells us that the polarization is correlated with the 
temperature perturbations. The correlation power Cf^ being the product of the two, exhibits oscillations 
at twice the acoustic frequency. As in the case of the damping, the predicting the precise value requires 
numerical work pO| ] since r changes so rapidly near recombination. Nonetheless the detailed predictions 
shown in Fig. || bear these qualitative features. 

The acoustic polarization and cross correlation has recently been detected by the DASI experiment p5[ . 
Like the damping scale, the acoustic polarization spectrum is uniquely predicted from the temperature 
spectrum once fifc/i^, i^mh^ are specified. Polarization thus represents a sharp test on the assumptions of the 
recombination physics and power law curvature fluctuations in the initial conditions used in interpreting the 
temperature peaks. 

These acoustic peaks in the polarization appear exclusively in the EE power spectrum due to the azimuthal 
symmetry of the plane wave fluctuations. During the break down of tight coupling that occurs at last 
scattering, any gravitational waves present will also imprint a local quadrupole anisotropy to the photons 
and hence a linear polarization to the CMB These contribute to the BB power and their detection 
would provide invaluable information on the origin of the fluctuations (see e.g. p^). Specifically, in simple 
inflationary models their amplitude gives the energy scale of inflation. The gravitational wave amplitude h 
oscillates and decays once inside the horizon, so the associated polarization source scales as h/r and so peaks 



at the / w 100 horizon scale and not the damping scale at recombination (see Fig. This provides a useful 
scale separation of the various polarization effects. 



5. SECONDARY ANISOTROPY 

5.1. Gravitational Secondaries 

As the photons travel from the last scattering surface to the observer, they traverse the large-scale structure 
of the universe leaving subtle imprints in the temperature and polarization fields. 

There are two main gravitational effects. Intervening mass along the line of sight gravitationally lenses the 
CMB photon trajectories |6|] and hence distorts both the temperature and polarization anisotropy fields. 
The photons are deflected according to the angular gradient of the potential projected along the line of sight, 

cj){h) = 2j d^i-^__i$(Dn,,7). (25) 

Because surface brightness is conserved in lensing, the deflection simply remaps the observed fields as 

X{n)^X{n + W(l)), (26) 

where X e {&, Q, U}. The typical deflection angle is of order a few arcminutes but the lines of sight are 
coherently deflected across scales of a few degrees. Since the coherence scale of the acoustic features is larger 
than the deflection angle, the lensing effect can be calculated by Taylor expanding Eqn. ( p6| ) (e.g. |6^). The 
result is a product of fields so that in harmonic space the modes are coupled to each other across a range 
Al « 60 set by the coherence of the deflection. Heuristically, lensing distorts the hot and cold spots formed 
by the acoustic oscillations and hence the mapping of fc to Z. 

In the temperature power spectrum, this mode coupling smooths the acoustic peaks slightly |70[ . For the 
polarization, the remapping not only smooths the EE-peaks but actually generates B-mode polarization [0 . 
Remapping by the lenses preserves the orientation of the polarization but warps its spatial distribution and 
hence does not preserve the symmetry of the original E-mode. Gravitational lensing represents a fundamental 
obstacle to detection of low amplitude gravitational waves from inflation. 

Because the lensed CMB distribution is not linear in the fluctuations, higher order statistics are promising 
probes of lensing effects ||7|, 73, |2[ In particular, the coupling of multipoles separated by Al < 60 can 



be used to construct a minimum variance estimator of the deflection potential (j}{n) out of pairs of moments 
iffsf . These quadratic estimators are close in performance to the optimal but more complicated maximum 
likelihood estimator |7^. Since small-scale polarization anisotropy is otherwise free of cosmological i?-modes, 
most of the signal-to-noise in the reconstruction lies in the pairing of iS-modes to neighboring lensing- 
generated i?-modes | |77| . This correlation can in principle be used to map the deflection potential down to 
the 10' scale. This mapping can in turn be used to measure the growth of structure from — 10"^ and 
hence the properties of the dark matter, including any component from massive neutrinos. It can also be used 
to remove much of the contamination to the gravitational wave _B-modes and bring the detection threshold 
for the energy scale of inflation down to a i^i « few 10^^ GeV pof . 

Gravitational potentials can also change the temperature anisotropy through gravitational redshifts. 
Density perturbations cease to grow once the dark energy dominates the expansion. As in the case of the 
matter-radiation transition, the gravitational potentials must then decay. The evolution of the gravitational 
potential under a smooth component of dark energy is given by 



|[l-u;e(a)]f7e(a)$ = 0. (27) 



d\na 2 



With the appropriate replacements, this equation also applies to smooth components of dark matter, e.g. a 
l^ht neutrino. A scalar field candidate for dark energy is smooth out to the horizon scale at any given time 

As we have seen, a decay in the gravitational potential causes an effective heating of the photons in a 
gravitational well. Like the matter-radiation analogue, it is intrinsically a large effect since the net change 



due to the decay is 5 times the Sachs- Wolfe effect of ^/3. However since the opposite effect occurs in 
underdense regions, the contributions are canceled as photons traverse many crests and troughs of the 
potential perturbation during the matter-dark energy transition. The effect, called the (late-time) integrated 
Sachs- Wolfe (ISW) effect Q , then appears only at large scales or low multipoles (see Fig. ||) . 

The ISW effect is the most direct signature of the dark energy available in the CMB and is very sensitive 
to its equation of state and clustering properties if it deviates far from a cosmological constant {we > —1/2, 
see Fig. unfortunately given the large-scale nature of the effect, the limited number of samples of the 
low order multipoles, or "cosmic variance" , prevents a precise determination from the power spectrum. A 
promising technique is to isolate the ISW effect through its cross correlation with reconstructed lensing maps 
from the CMB @|8|, |7|l. 

5.2. Scattering Secondaries 

The universe is observed to be ionized out to 2 « 6 and is thought to have undergone reionization sometime 
between between 7 < z < 30 (see Q for a review). Consequently a minimum of a few percent of the photons 
have rescattered since recombination (see Eqn. o). The main effect from rescattering of the photons during 
reionization is a uniform suppression of the peaks by e~'^ as anisotropy is destroyed by the randomization of 
directions in scattering. Since this can be confused with a change in the initial amplitude of fluctuations, it is 
important to resolve this ambiguity for the study of initial fluctuations and the growth rate of fluctuations, 
which itself probes the properties of the dark energy and dark matter through Eqn. (|2^). Reionization is 
interesting in its own right since it can us about the sources of ionizing radiation from the first astrophysical 
objects formed in the universe. 

The most promising means of isolating the reionization effect through the CMB is from the generation of 
large scale polarization during recombination. The rescattered radiation becomes polarized since temperature 
inhomogeneities become anisotropics by projection (see Eqn. |8|), passing through quadrupole anisotropy when 
the perturbations are on the horizon scale at any given time (fcZ3* k, krj k,1 = 2). The result is a peak in 
the iJ-mode power spectrum shown in Fig. |l|. I n a perfect, foreground-free world, even the minimal signal 
is within reach of the MAP [|4| and Planck satellites and can be used to isolate the reionization epoch 

The dissipation of the acoustic oscillations at small scales also provides a window through which to see 
subtler effects. The main contribution beyond the damping tail comes from the rescattering of the CMB off 
of hot electrons in clusters of galaxies, called the Sunyaev-Zeldovich effect [Q. Inverse Compton scattering 
represents a net transfer of energy between the hot electron gas and the cooler CMB. It leaves a spectral 
distortion in the CMB where photons on the Rayleigh- Jeans side are transferred to the Wien tail. The 
effect in individual clusters is now routinely measured and is an important means of studying the physics 
of clusters and the dark energy . The effect of unresolved high redshift clusters may have already been 
detected in the CMB by the CBI and BIMA experiments, albeit at an amplitude that is somewhat higher 
than expected [ p5| |9l[ |. Confirmation of the nature of this small scale excess will require measurements at 
multiple frequencies to test its inverse Compton nature. 

There are a host of smaller effects due to Doppler shifts of the moving ionized gas (see e.g. Q). These 
present a significant challenge to detect and interpret and are beyond the scope of this review. 



6. DISCUSSION 

The tremendous experimental progress in CMB anisotropy studies over the last few years has helped establish 
a standard paradigm for cosmology at intermediate epochs but has simultaneously raised questions about 
the physics at the two opposite ends of time. Simple inflationary models of the early universe have so far 
passed the test of the acoustic peaks. They predict the near scale-invariant initial curvature fluctuations in a 
spatially flat background that have now been observed. Moreover, spatial flatness appears to be maintained 
today by a mysterious form of missing or dark energy. 

The first step into the future will come with the release of data from the MAP satellite in early 2003 which 
in the full course of the mission will provide the definitive measurement of temperature anisotropy across 



the whole sky down to a quarter of a degree. The MAP satelhte wiU test the spectrum of initial conditions 
beyond the simple constraints on the power law index which can be extracted today (e.g. |9^). These 
advances will come with the increased precision in the temperature measurements and determination of the 
acoustic peaks in the polarization and cross correlation. Measurement of the polarization peaks will not only 
eventually double the amount of statistical information that can be extracted from the CMB but can also 
provide a sharp distinction between small changes in the dynamics of the plasma, which affect temperature 
and polarization differently, and small deviations from scale invariance of in the initial conditions (e.g. [Q), 
which affect the two alike. Studies of the Gaussianity of the acoustic fluctuations will also provide a strong test 
of the inflationary model. The increased precision on the matter-radiation ratio will sharpen considerably 
the constraints on the dark energy. MAP should also be able to detect or place limits on the amount of 
rcionization in the universe through a large angle bump in the polarization and hence remove a central 
ambiguity for dark energy and dark matter studies from the peaks. 

The next generation of experiments dedicated to fine-scale, secondary structure in the anisotropy will 
enable new tests of the dark energy from structures in the universe such as galaxy clusters. They will also 
provide a new handle on rcionization through the Doppler shift of CMB photons off of moving structures. 

Further down the road lies the milestone of the gravitational lensing of the peaks. Detection will not 
only require more sensitive instruments and higher angular resolution but also subtraction of galactic 



and extragalactic foreground contaminants |96[ |97|. In principle higher order statistics can be used 
to reconstruct the potential field in projection with a large gains in sensitivity available if the B-modes 
from lensing can not only be detected but accurately mapped fn\ . From this reconstruction the dark matter 
and dark energy dependent growth of structure can be measured. For example gravitational lensing can 
enable tests of the scalar-field hypothesis for the nature of the dark energy through cross-correlation with 
the integrated Sachs- Wolfe effect [fTSj . 

The ultimate future milestone for the CMB is the detection of gravitational waves from inflation through 
_B-mode polarization. Detection would represent strong evidence for the inflationary model and pin down its 
energy scale. The window of detectability can be extended to energy scales above a few 10^^ GeV but only 
if the i3-modes from gravitational lensing are removed from a direct reconstruction. 

These future milestones will require much experimental effort to achieve and much theoretical effort to 
interpret. Nevertheless they provide the hope that the next decade of studies of CMB temperature and 
polarization anisotropy will be as fruitful as the last. 
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